The spectrum of particles accelerated in relativistic, collisionless shocks 
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We analytically study diffusive particle acceleration in relativistic, collisionless shocks. We find a 
simple relation between the spectral index s and the anisotropy of the momentum distribution along 
the shock front. Based on this relation, we obtain s = (3/3„ — 2(3u0d + Pd)/{Pii, ~ Pd) for isotropic 
diffusion, where /3u (/3d) is the upstream (downstream) fluid velocity normalized to the speed of light. 
This result is in agreement with previous numerical determinations of s for all {l3u,Pd), and yields 
s = 38/9 in the ultra-relativistic limit. The spectrum-anisotropy connection is useful for testing 
numerical studies and for constraining non-isotropic diffusion results. It implies that the spectrum 
is highly sensitive to the form of the diffusion function for particles travelling along the shock front. 



Diffusive (Fermi) acceleration of charged particles in 
collisionless shocks is believed to be the mechanism re- 
sponsible for the production of non-thermal distributions 
of high energy particles in many astrophysical systems 
. This process is believed to play an important role in, 
for example, planetary bow shocks within the solar wind, 
supernovae remnant shocks driven into the inter-stellar 
medium Q, jets of radio galaxies 0], gamma-ray bursts 
(GRB's) and possibly shocks involved in the forma- 
tion of the large scale structure of the miiverse j^]. This 
phenomenon is common to shocks with widely differing 
physical characteristics (e.g. velocity and length scale), 
as evident from the above examples. 

The Fermi acceleration process in shocks is still not 
understood from first principles (see, e.g., (3] for a dis- 
cussion of alternative shock acceleration processes) . Par- 
ticle scattering in collisionless shocks is due to electro- 
magnetic waves. No present analysis self-consistently cal- 
culates the generation of these waves, the scattering of 
particles and their acceleration. Most analyses consider, 
instead, the evolution of the particle distribution adopt- 
ing some Ansatz for the particle scattering mechanism 
(e.g. diffusion in pitch angle), and the "test particle" ap- 
proximation, where modifications of shock properties due 
to the high energy particle distribution are neglected. 

This phenomenological approach proved successful in 
accounting for non-thermal particle distributions inferred 
from observations. The theory of diffusive particle accel- 
eration in non- relativistic shocks was first developed in 
1977 [1 nil 0. Diffusive (Fermi) acceleration of test 
particles in non-relativistic shocks was shown, in partic- 
ular, to lead to a power-law distribution of particle mo- 
menta, (Pn/dp^ ocp^*, with 

s = SPu/iPu - Pd) ■ (1) 

Here, /3„ {j3d) is the upstream (downstream) fluid velocity 
normalized to the speed of light. For strong shocks in an 
ideal gas of adiabatic index F ~ 5/3, this implies s — A 
(i.e. p^d^n/dp^ ocp^^), in agreement with observations. 

Observations of GRB afterglows lead to the conclusion 
that the highly relativistic collisionless shocks involved 
produce a power-law distribution of high energy parti- 
cles with s = 4.2 ± 0.2 ^ld\. This triggered a numerical 



investigation of particle acceleration in highly relativistic 
shocks JTj. The values of s were calculated, in the "test 
particle" approximation and assuming velocity angle dif- 
fusion, for a wide range of Lorentz factors and several 
equations of state [TTlll2 . IT^ and the references therein]. 
In particular, s was shown to approach the value 4.2 for 
large Lorentz factors, in agreement with GRB observa- 
tions. The study of particle acceleration in relativistic 
shocks is of interest to many other high-energy astro- 
physical systems as well, e.g. jets in active galactic nu- 
clei 01 and in X-ray binaries (micro-quasars) [H, and 
may be relevant for the production of ultra-high energy 
cosmic-rays 16]. 

The analysis of shock acceleration is more complicated 
in the relativistic case than in the non-relativistic case, 
mainly because the particle distribution near the shock 
is highly anisotropic. Due to this difficulty, only approx- 
imate numerical analyses (using Mote-Carlo simulations 
or eigenfunction methods) are available for relativistic 
shocks. In particular, an analytic expression for s ex- 
tending Eq. JQ) to the relativistic regime is unavailable. 
In this paper we present an analytic study of diffusive 
particle acceleration in relativistic shocks, under the test 
particle and velocity angle diffusion approximations. 

/. Formalism. Consider a shock front perpendicular to 
the z-axis, where the fluid flows in the positive z direction 
upstream (z < 0) and downstream (z > 0). For diffusion 
in the direction p of rest frame momentum p, the station- 
ary transport equation for particles with Lorenz factors 
well above the shock Lorenz factor, is 17] 

7i(A + fJ-i)dzMfJ.i,Pi,z) = d^ji{fii,pi,z) , (2) 

where upstream/downstream indices i G {u,d} will be 
written only when necessary, 7 = (1-/3^)^^/^ is the 
Lorenz factor, and f{fj,,p, z) is the (Lorenz invariant) par- 
ticle density in the phase space of = cos(p ■ z), p, and 
the distance from the shock front z as measured in the 
shock frame. The flux in momentum space is 

j{fi,p,z) = D^f,{fi,p,z)d^f{fi,p,z) , (3) 

where i?^^ is the diffusion function, and oc {1 — ^P) 
for isotropic diffusion. Assuming that i?^^ is separable 



2 



in the form Di{ii)D2ip, z), Eq. Q may be written 

iP + ^,)^rfi^l,p,T) ^ d^[{i - ^l^)D{^l)^^,f{^l,p,T)] , (4) 

where -D(^) = Di{^)/{1 — and t = z D2{p, z)/'-f. 

Next, wc incorporate boundary conditions. Continuity 
across the shock front imphes 



U ; Pu 



0) = fd{pd,Pd,Td = 0) 



(5) 



where upstream and downstream quantities are related 
by a Lorenz boost of velocity Pr = {Pu — Pd)/{^ — PuPd)', 
Pd = 7rPu(l+/3rMu) and = +/3r)/(l +/3rA*«)- Par- 
ticle injection near the shock, and diffusion of particles 
into the far downstream only, imply that 



lim /„ = and lim = 



(6) 



where /oo is a constant. 

The spectral index s was previously calculated, by nu- 
merically finding some eigenfunctions of Eq. that sat- 
isfy the boundary conditions of Eq. ©, and approxi- 
mately matching the upstream and downstream solutions 
El III m, or by Monte Carlo simulations [Hill. Ref. 
[l3| has also exploited the relation 



s = 3 - 



HPret) 

ln{Ef/E.,) 



(7) 



where (Pret) is the (flux averaged) probability of a par- 
ticle crossing the shock downstream to return upstream, 
and {Ef/Ei) is the average energy gain per cycle. In the 
ultra-rclativistic limit, Pu = 1 and Pd = 1/3, such studies 
have converg ed on the value s = 4.22 ± 0.02 [Hill El. 

//. Anisotropy- Spectrum Connection. The particle 
drift downstream implies that the distribution near 
the shock front is anisotropic^ more particles moving 
downstream than upstream. The distribution is nearly 
isotropic in the non-relativistic limit, where Eq. is 
approximately recovered by demanding that fd{Td = 0) 
is isotropic jig- The lack of a characteristic momentum 
implies that the spectrum remains a power-law in the 
relativistic verified numerically [ill Il3l| . 

In a steady state, the particle distribution / is station- 
ary parallel to the shock front, implying that d^j = for 
fx = —P, as evident from Eq. Lead by this obser- 
vation, we expand f{^,p,T) and -D(/Lt) around /z = —P 
(assuming smooth functions in this region, see e.g. [I2j|): 



fi{pi,Pi,Ti) = a''^\n) + a\^'{Ti){fi^ + Pi) 



(8) 



+ p,y + ■ ■ 

D^ifi^) = + d['^ ifi, + P,) -f + p,f + ■■■ 

Using Eq. |SJ|, one can relate the coefficients aj{T — 0) 
on both sides of the shock front. For the two lowest order 
terms, these relations yield an explicit expression for s; 



rd 



Pu - Pd ' 



(9) 



where we have defined 
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] 



DD2f 



(10) 



T=0,fi=-(3 



The physical significance of ai and oq may be demon- 
strated as follows. For rest frame velocity angle diffu- 
sion, the mass and momentum flux are generally not con- 
served, and particle energy is conserved only in the rest 
frame. Conservation of particle number and rest frame 
energy yield the one-dimensional continuity equation. 



df J d^if{^i,p,z)+d^' J dnfj,f{n,p,z) ^0 , (11) 

where t' and z' are rest frame variables. Transforming to 
shock frame variables, we find that in a steady state 



dfi{P + fj.)f{fj.,p,z) = Cp' 



(12) 



where C is a constant determined by the boundary con- 
ditions; Cu = upstream and Cd = 2Pdfoo > down- 
stream 19j. If we define the convection towards the up- 
stream and towards the downstream as 



G+ — I dug and G- — 



dfii-g), (13) 



then and G_ are non-negative, G-^- — G_ upstream, 
and G+ = G_ + Cd downstream. In addition, 



drG+ = drC- = -ai{T)D{fi = -p)j- 



(14) 



reflecting the fact that as particles turn around from 
heading upstream to heading downstream or vice versa, 
they must diffuse through a state where they propagate 
parallel to the shock front. Finally, oq = /(/^ = — /?) is 
the normalization of / along the shock front. 

Next, we exploit the stationary particle distribution 
parallel to the shock front. Substituting Eq. © into 
Eq. Q for fi = — /?, we find a simple connection between 
the coefficients ai and 02, valid for any r: 



a2(T) = -aiir) Pj^ + 



1 di 
2d^ 



(15) 



When extrapolated to t = on both sides of the shock 
front, this result may be combined with the relation be- 
tween 02"^ and 02^^^ (cf. Eq. and Eq. [H]) and with 
Eq. , to yield a relation between r„ and 

Tuiru + du + Pu + Pd) = rd{rd + dd + Pu + Pd) , (16) 

where di = T^^'^i'V^*'* is a measure of the deviation 
from isotropic diffiision, for particles moving almost par- 
allel to the shock front. Alternatively, Eq. jnj and 
Eq. 1)16(1 may be combined to yield expressions for s as a 
function of on only one side of the shock. 



sird) 



2{Pu - Pd) 



[Pu+Pd + 2rd + du± 



(17) 



Vif^u + Pd + 2rd + duY + 4rd(dd - d„) 
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and 

siru) 



[-{/3u+f3d + 2r„ + dd) ± (18) 



2{Pu - Pd) 
ViPu +Pd + 2r„ + dd)^ + 4r„(d„ - dd) 



For isotropic diffusion, where Di{^) is constant, di = 0, 
Eq. lfTB|) simplifies to + Vd + f3u + Pd = 0, and 



2rd + l3u + l3d {-2ru) - f3u - f3d 



I3u-Pd 



Pu-Pd 



(19) 



Tlie above analysis can be repeated by expanding the 
distribution function in any other frame. For example, 
expansion around « in the shock frame yields 



sirs) 



1 



Pu+Pd 



duPu - ddPd ± 



(20) 




Pu+Pd 



+ duPu-ddPd] +2rs{dd-du) 



where di = di/(2/3^ — 2/3^), and Ts = a^^'^ /a\^'' is related 
to r„ and through = {PuTd - Pdru)/{Pu - Pd)- For 
isotropic diffusion, Eq. (|2()() simplifies to 



l + 2rJ{P^ + Pa) . 



(21) 



///. Analytic expression for s. The boundary condi- 
tions at r ^ ±oo were not used in the analysis leading 
to Eqs. (|17I21|I . but are needed in order to determine the 
anisotropy parameters ai/oo. In general, the convection 
gradually decreases at larger distances from the shock 
front. Eq. ^ thus implies that a^"' > and aj''^ < 0. 
For isotropic diffusion, this yields the constraint 



S^so > {Pu+ Pd)/{Pu- Pd) 



(22) 



More important consequences of the above connec- 
tion derive from cases where s is known; (i) the non- 
relativistic limit of Eq. J^l; and (ii) the limit of infinite 
compressibility Pd = 0, where (Pret) = 1 and thus 

siPu > 0, /3d = 0) = 3 . (23) 

For isotropic diffusion, Eqs. H19|) and (|23|l imply 

af^/a'if^ = Pu-^Pd + PdAdiPu^Pd) and (24) 



(u) , (m) 



2Pu + ^Pd + PdAuiPu,Pd) 



The unknown functions Ai are at least second order, even 
pof functions of /3„ and Pd- Hence, in the non-relativistic 

limit, a[''^/alf^ ~ Pu-Pd/2 and a^'Voo"^ - -2Pu-Pd/2. 

In the relativistic case, the shock front distribution is 
far more isotropic in the downstream frame than it is in 
the upstream or the shock frames. In particular, ai/aQ 
includes (at least) third order terms (in /?„ or Pd) up- 
stream, but could be of first order downstream. Using 



only the first order terms of oi/ao downstream (i.e. as- 
suming Ad = 0) yields a simple expression for s. 



{3Pu-2PuPj+P!i)/{Pu-pd) 



(25) 



which is in excellent agreement with numerical studies 
llSs JJJ over the entire relevant range of /3„ and Pd- This 
is demonstrated in Figure for three different equations 
of state. In the ultra- relativistic limit, Eq. H25|l implies 



,(/3„ ^l,Pd-^ 1/3) = 38/9 = 4.222 . . 



(26) 



in excellent agreement with Refs. The ex- 

pression a'^^ I a^Q^ = Pu—Pdl2 itself agrees with the down- 
stream distribution calculated numerically 12, Figure 3] 
even in the ultra-relativistic case, although we can not 
prove that this first order expansion is sufficient. 



FIG. 1: Comparison between Eq. I|25|l (curves) and previ- 
ous numerical results []3 (symbols). The spectral index s is 
shown as a function of /3u7ti for three different types of shocks; 
a strong shock with the Jiittner/Synge equation of state (solid 
curve vs. crosses), a strong shock with fixed adiabatic index 
r = 4/3 (dashed curve vs. x-marks), and a shock in a rela- 
tivistic gas where /3„/3d = 1/3 (dash-dotted curve vs. circles). 
For a description of these equations of state, see Ref. pl|. 

lY. Other consequences- Eqs. (|17I21|) directly relate 
the spectral index to the distribution of particles trav- 
elling nearly along the shock front. Note that Eqs. 
((T7|l . ifT^ and (|2nil hold for any form of the diffusion 
function. These results provide an important consis- 
tency check for numerical methods such as used in Refs. 

[afiiiiiEiii. 

In some cases, oi/ao can be read directly from the 
calculated particle distribution. For example. Figure 3 
of Ref. [l3l shows the shock front distribution in both 
downstream and shock frames, for isotropic diffusion in 
an ultra-relativistic shock with (/3„ = 0.995, Pd — 0.328). 



The figure yields a^^'/a\^' ~ 0.85 and a^f /a[''^ 
corresponding to s ~ 4.2 according to both Eq. H19|l and 
Eq. H21|) . For the same shock and non-isotropic diffusion 



2.1, 
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with Dili) oc ifi'^ + 0.01)-i/2, Figure 5b of Ref. [13 
implies that a[^^ /a'jj '' ~ 1.9. Using Eq. if^ . this also 
corresponds to s ~ 4.2. 

The spectrum-anisotropy connection may be used to 
associate each eigenfunction k of Eq. with its cor- 
responding spectral index Sk- For isotropic diffusion 
in the ultra-relativistic limit, the shock front distribu- 
tion coincides up to a 10% accuracy with the first up- 
stream eigenfunction /i where in the shock frame 
/i oc (1 - PufJ^s)'" exp[-(l + /^s)/(l - (3u^J.s)]■ According 
to Eq. this gives si = (2 + /3„ - Pd)/{Pu - /3d) = 4 
in the ultra-relativistic limit. 

For non-isotropic diffusion, Eqs. I|17|l and H23() imply 
that, to lowest order in f3d, 



f3u - f3d/2 - du/2 
Pu + [dd - du)/6 



(27) 



For the non-isotropic diffusion function solved numeri- 
cally 12], D{p) oc (/i^ -I- 0.01)~^/^, the approximation 
Eq. (|27|l is in poor agreement with the numerical data, 
although in the ultra-relativistic limit it yields s = 4.26, 
similar to the result s = 4.21 of Ref. "T5|. 

It is interesting to note that Eq. (|25|l may be written 
in the form s^so = 3 - (-3/3^ -I- 2PuPj - PD/iflu - Pd). 
Comparing this result with Eq. Q may suggest that 



(P.et) =exp(-3/3d + 2/3„/3j-/3;^) 
and {Ef/E^) =ei^-p{l3u- 13d) 



(28) 



Although other choices that conserve In(Pret)/ ^'^{Ef /Ei) 
are possible, Eq. H28|l does agree well with numerical 
studies of the ultra-relativistic limit, where (Pret) = 
0.439 ± 0.007 and (Ef/E^) = 1.97 ± 0.02 

V. Summary. We have analytically studied the spec- 
trum of test particles accelerated by an arbitrary rela- 
tivistic shock in the diffusion limit. A simple relation 
was shown [Eqs. p7|l - 12H) ] to exist between the spectral 
index s and a measure of the particle anisotropy, ai/a^. 
For isotropic diffusion, the lowest order terms in ai/oo 
downstream yield an expression for s [Eq. H19|l ]. which 
is in excellent agreement with previous numerical studies 
(Figure over the entire relevant range of and (3d- In 
the ultra-relativistic limit, it yields s = 38/9, in agree- 
ment with previous studies and with GRB observations. 

The spectrum-anisotropy connection in Eqs. H17|) . H18|) 
and H2U|) holds for any diffusion function D{fi). It indi- 
cates that s is particularly sensitive to the form (first 
derivative) of D{fi) for upstream and downstream par- 
ticles travelling along the shock front. This connection 
is also independent of the test particle approximation, 
providing a useful tool or at least a consistency check for 
various studies in the diffusion limit, including a future 
self-consistent calculation of particle scattering, acceler- 
ation, and electro-magnetic wave generation. 
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